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Atatmci 

Underwater towing experiments wore carried out 
with a rectangular airfoil of aspect ratio 5,3 at 4° 
and 0° angles of attack and at chord-based Reynolds 
numbers between 2nlO s and 7, 5x10 s , ouantitative mea- 
surements by means of tho hydrogen bubble technique 
indicated lower peak swirl velocities in the range 
of 100 to 1000 ciiord lengths downstream than have 
been measured in wind tunnel or flight tests, The 
maximum circumferential velocity decayed as 
whereas the turbulent eddy viscosity, n.j., increased 
as t 3 / 4 . This behavior and other known rates of vor- 
tex decay are explained in teims of an analytical 
solution for the vortex problem with timo varying 
eddyviscosity.lt is shown that this case corresponds 
to nonequilibrium turbulent vortex flow which has 
other self similar solutions, such that the rate of 
vortex docay is greater than fl/2 when the eddy vis- 
cosity increases with time, Conversely, the decay 
proceeds at a rate loss than t -1 ' 2 when decreases 
with time, 

flotation 

a = core radius 

b « wing span 

c B wing chord 

AH ■* head loss in boundary luye* 

m “ exponent for time vurying eddy viscosity 

p « pressure 

a fluid pressure m r 4 “ 
p(0) “ pressure on the yortex axis 

p ■ rate of turbulent energy production 
r ■ radius measured from vortex axis 
s “ tangential shear deformation, s * 3v/9r - v/r 
t “ time 

At “time difference in velocity measurement 

T ■ time parameter, T ■ r^t/c 2 

u ■ radial velocity component 

Uoo ■ towing speed 

v ■ tangential velocity component 

v a « tangential velocity at r « a 

-jjT^r » Reynolds stress 

w = ajcial velocity component 

w 0 ■ axial velocity on vortex axis 
t ■ downstream distance behind airfoil 

a * angle of attack 

0 “ constant of proportionality, Eq. (10) 

P « circulation 

«■ circulation at core radius 
fuo » circulation of wing 
0,AO a angle in velocity measurement 
v ■ kinematic viscosity 

F c «r effective viscosity 

P-p « eddy viscosity 

(Pf/a) a B eddy viscosity at core radius 
(Pf/p) o* eddy viscosity level after roll up 
P » density 

t * turbulent shear stress 
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introduction 

With the advance of the heavy transport jet 
planes such rt tho 747, and more recently tho DC-10 
and L-1011, a rather intense theoretical and experi- 
mental effort lias boon directed toward the problem 
of ’'wake turbulence," This so-called turbulence, also 
known as turbulent wako, consists of the remains of 
a system of counter rotating wing tip vortices which 
are shed from the aircraft as part of the generation 
of lift, When a following plane enters the vortex 
wake it can, depending upon its site, be rolled out 
of control or, in the more severe cases, even suffer 
serious structural damage, 

The intensity and longevity of vortex wakes in- 
crease with tho weight of the aircraft generating 
tho wake, but in general, the severity of the upset 
is proportional to tho swirl momentum left in an ef- 
fective cross sectional area of the wake at tho timo 
of encounter, Eince the wake vortices dissipate in 
time as a result of turbulent and atmospheric effects, 
an encounter with an "old" wake is expected to in- 
duce less severe roll on the following craft than one 
with a newly formed wake, Present aviation rules call 
for several miles separation between planes on ap- 
proach, so that possible encounters would occur in 
tho far field of the wake where the separation be- 
tween planes is on the order of 10 3 or more mean 
chord lengths, There are, however, strong incentives 
toward closer aircraft spacing and hence toward in- 
creased hazards from encounters with vortex wakes. 

In several recent wind tunnel studies, most 
notably those of Chigier and Corsiglia, Jo 2 Logan, 3 
and Mason and Marchman, 4 tangential and axial velo- 
city components of tho vortex field wero measured 
after roll up. Subsequent velocity' distributions and 
initial rates of decay of the flow field were obtain- 
ed for downstream distances up to z/c » 30, and can 
thus be regarded as data from the near and interme- 
diate fields of the wakes, 

Although flight measurements of maximum tangen- 
tial vortex velocities acquired as far back as 7000 
chord lengths have been reported by McCormick 
ot al., 5 no such data are available for comparison 
from model tests where the dynamics of vortex decay 
in tho far field can be studied under well controlled 
conditions. Since in-flight encounters with vortex 
wakes are most likely to occur in the far field, a 
better understanding of the mechanics and the struc- 
ture of wake vortex flow in the late stage is re- 
quired, 

The purpose of the present paper is to report 
on a series of underwater towing experiments which 
allowed measurements in the far wake region of a 
rectangular NACA-0015 airfoil of 0.915 meter span 
with rounded wing tips and aspect ratio 5,33. Chord- 
based airfoil Reynolds numbers ranged from 2,2x10 s 
to 7,5x10 s , and angles of attack were set at 4° and 
B°. Use of the hydrogen bubble technique** permitted 
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*This technique uses hydrogen bubbles generated by 
electrolysis on a metsl wire to mark the flow. 



visual observations! of this development of the turbu- 
lence structure in the core and the potential regions 
during vortex decay, Simultaneously, the flow was ’ 
photographed on 16 raa movie film from an axial and a 
crosswise viewing direction. Quantitative measure- 
ments were made from these films of the tangential 
velocity profiles and maximum axial velocity, Tlius, 
rates of vortex grout!) measured at the core and de- 
cay of the maximum tangential velocity were estab- 
lished as a function of chord-based downstream dis- 
tance, 

In the last section of this paper, a ciosed- 
forts solution is obtained for the vortex decay with 
time varying eddy viscosity, which affords determi- 
nation from the data of eddy viscosity levels in the 
vortex flow at given points in time. This information 
is then used to compute radial distributions of Cor- 
fu lent shear stress and rates of production of cur- 
Iiulence energy, 


The towing tests wore conducted in the Lockheed 
Underwater Missile facility (4,6 n x 4,6 ta x 5S m) 
with the airfoil mounted to a towing carriugo by a 
streamlined strut (see Fig. 1), The hydrogen bubble 
generating wire consisted of three twisted platinum 
wires, 4.6 h dia, and was mounted within the tank at 
an angle to the vertical and bolcw the right wing 
tip. A S kw light source mounted above the tank was 
capable of illuminating the plane of the platinum 
wire witn an intense light sheet, 1,2 m x 0,3 tn in 
cross section. Two cameras were mounted at the ele- 
vation of the airfoil, one of which was lined up 
within the plane of the platinum wire; the view of 
the other camera was pointed in the axial direction 
by use of a mirror, 

Side view and axial view film recordings of the 
right vortex were made in two modes, one in which 
hydrogen bubbles were generated at the rato of 3/sec 
from the moment of passage of the airfoil until the 
vortex* flow field was well marked. This mode of op- 
eration, illustrated in Fig, 2, yielded excellent 
photographs for further visual study. I 1 ' the second 
mode, lines of hydrogen bubbles were only pulsed 
while the vortex core crossed the wire, so that they 
were not obscured by previously marked fluid, These 
lines, also visible in Fig 2, served for measure, mont 
of the tangential velocity profiles, For the purpose 
of velocity measurements, the ago of the vortex was 
determined by the time elapsed between passage of 
the airfoil and the instant at which tho vortex core 
crossed the wire. The center of the core, or vortex 
filament, was usually well mark)*! by hydrogen bubbles 
left in the water from the previous run. 


Tangential Velocity and Circulation 

The method of obtaining velocity profiles by 
means of the hydrogen bubble technique in a linear 
flow is described in detail by Schraub et al,° In 
the present work, radial velocity distributions were 
deduced from film by measuring f~om the center of 
the vortex the polar coordinates of the hydrogen 
bubble time line next to the wire, whereby the down- 
ward motion of the vortex system, given by r<j2nb, 
was neglected. Since the time delay, At, between 
time lines is known from the film frame speed, tho 
tangential velocity variation is given by 


v(r) u CD 

This procedure is indicated schematically in Fig 3. 
After computation of v(r) for each side of tho vor- 
tex, the results were averaged in order to represent 
tho vortex by one velocity profile, 

Typical velocity profiles, normalised on the 
towing speed, U^, are shown in Fig 4 for 4° and B e 
angles of attack. The radius is normalised on the 
wing span, b. The velocity profiles decay with down- 
stream distance in a manner quite similar to a lami- 
nar vortex! this point will be discussed further in 
a later section, The reduction in maximum tangential 
velocity proceeds with a simultaneous enlargement of 
the core, outside the core, the tungentiul velocity 
approaches the envelope of the potential vortex, but 
it should be noted that due to viscous or turbulent 
shear forces, the flow field decays everywhere simul- 
taneously. 

Fig S shows the decay of the maximum tangential 
veloci* • for the range of downstream d. stance inves- 
tigated in this study, For comparison, we_also list 
the flight tost data of McCormick et al,, which was 
normalised appropriately, The effect of angle of at- 
tack upon v u /U , shown in the flight test data* and 
also noted in F.efs, 2, S, and 7 for 'short distances 
behind the wing, has disappeared from our data for 
the downstream distances investigated, As will be 
shown in the discussion of turbulent contributions 
to vortex decay, higher levels of turbulence, or 
eddy viscosity, were produced at higher angles of 
attaek. The result was an increased rate p£ decay of 
v a , which tended to diminish the differences in maxi- 
mum tangential velocities caused by tho dependence 
upon ot, The data indicate, furthermore, that for t/c 
above 300, the rate of decay of the maximum swirl 
velocity approached a slope of -7/3. This rate is 
considerably higher tlrnn the t "4/4 dependence exhi- 
bited by the flight test data of McCormick et a,, 5 
for large distances behind the airfoil. Using the 
notion that turbulent diffusion is proportional to 
the root mean square of the turbulent fluctuations, 
Brown 8 concluded that the turbulent vortex decay 
should follow as t" 1 ' , hence (t/c) -1 '' 3 , indicating 
a slower rate of vortex decay than that usually oo- 
served. The computations by Baldwin et al. 9 on tur- 
bulent vortices, however, showed that turbulent dif- 
fusion contributes significantly onj.y in regions of 
low shear. Brown's analysis, applied to the region 
of the core radius, does not adequately explain the 
vortex spreading rato, since the turbulent shear 
stresses are maximum just outside this region. Ex- 
amination of the present velocity data with respect 
to eddy viscosity levels within the trailing vortex 
showed that the high rates of decay and transport of 
core angular momentum outward were, in fact, due to 
the increase of the eddy viscosity with time. This 
discussion is deferred until a later section. 

The growth of tho vortex core region with down- 
stream distance is illustrated in Fig, ft. The indi- 
cated growth rate proportional to (c/c)*’ 1 ' 3 for the 
data at 4" agrees only coincidentally with the vor- 
tex spreading rate proposed by Brown, As mentioned 
above, the rate of radial transfer of angular mo- 
mentum is related to changes of the eddy viscosity 
with time. The higher growth rate observed at 8° 

•It is estimated that within the flight data of 
McCormick et al,, 8 the angle of attack varied by a 
factor of two. 
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angle of attack Is found to be the inverse of the 
rate of vortex decay deduced from Fig, 5, hit this 

condition need only be satisfied if we require that 
the product v(a,t)a(t) remain invariant with time. 
The faster spreading of the core angular momentum 
can be associated with higher turbulence levels in 
the vortex after roll up and with increased rates of 
turbulence production by ^hc mean shearing deforma- 
tion outstde tho core, 

Typical circulation profiles as a function of 
downstream distance are shown in Fig, 7. The trend 
toward self similarity, also evident in the velocity 
profiles, is well demonstrated by the data, over- 
shoot of the circulation was observed in some cases 
and may have been related to radial variations in 
tho rote of turbulent momentum transport, as dis- 
cussed by Govindruju and Saffman, 

Hoffman and Joubert 11 proposed a universal cor- 
relation of F(r)/P a vs log (r/u) . Fig, 8 shows the 
present circulation data in these coordinates, It is 
noted hero that this curve expresses merely the self 
similar property of the circulation profilos, but it 
cannot be regarded as a correlation of the flow on 
characteristic turbulence quantities, In the area of 
the core, tho centrifugal force field exerts a pro- 
nounced dumping effect upon the turbulent mixing 
process, thereby decreasing the length scale of tur- 
bulent momentum exchange. Donaldson* 2 pointed out 
that the use of eddy transport models docs not lend 
sufficient generality to the development of turbu- 
lence in aircraft trailing vortices. This short- 
coming accounts for the considerably narrowed extent 
of the loguritlmic region (0,6 < a < 1,3) when com- 
pared with the logarithiiic law of the turbulent 
boundary layer, 


associated with runs in which \>jiv increased in tho 
radial direction at the location of the core radius. 
For the remaining cases it was found that » r (v did 
not depond on r, 


Axial Velocity 

The importance of the axial core velocity has 
been previously discussed by Batchelor , 15 Brown, 8 
and Saffman, 16 As a result of these studies, it has 
become evident that the axial velocity in the core 
can be directed either toward tho wing or away from 
it, depending upon a balance between the viscous 
airfoil drag, which causes a wake to form, and tho 
swirl induced radial pressure gradient, which tends 
to propel the coro fluid as a Jet down the cortex 
core. Batchelor's formulation for the radial varia- 
tion of the axial velocity, expressed in coordinates 
of a stationary observer, is given by 

/ “» 

where 4(1 is the head loss due to the viscous boundary 
layer on the wing, If the viscous drag term is larger 
than tho integral term, as is the case at moderate 
anglos of attack, one observes a velocity defect in 
the core. However, the data of Chigier and Corsiglia 4 
show a reversal from a wake flow "in the core to an 
axial velocity excess, measured at a downstream dis- 
tance of z/c •< 9, when the angle of attack increases 
from 8° to 12°, Hence, at large angles of attack, 
the swirl induced pressure drop along the vortex 
axis immediately behind tho wing is strong enough to 
accelerate fluid from the wing tip area in the down- 
stream direction, despite the momentum loss incurred 
by passing through the viscous boundary layer, 


The value of the circulation at the location of 
maximum tangential velocity has been the point of 
much previous discussion. When the vortex decays 
under the effect of constant molecular viscosity, 
Lamb's 1 '- solution yields tho value r a /r M " 0,715, 

In the work of other investigators (see S?ef« 10), 
the value of P^/P,*, ranged from 0.37 to 0.6, and the 
model of ewindroju and Saffman 10 allows for 
r a /r M “ 1,2, In the recent vortex measurements with 
a laser velocimeter by Qrloff and Grant 2 at z/c ■ 2 
behind a rectangular wing, the coro circulations 
varied between 0.2 and 0,8 x H*, , depending upon tho 
angle of attack. The wind tunnel measurements of 
Corsiglia et al, 1 *’ resulted in r^/PcM Rs 0,4 - 0.5, 
Donaldson's 12 calculations indicate a slight de- 
crease with time of the initial circulation ut the 
core, whereas McCormick et al. proposed that when a 
decaying vortex maintains geometric similarity, 
r /r M » constant. This conclusion seems to be con- 
firmed by the data presented herein. Fig, 9 shows 
that tho circulation at the core radius remained 
essentially constant with downstream distunco at an 
average value of r a /Foo “ 0,67, The values plotted in 
Fig. 9 were determined by two different methods; 
with the first method, P_/Pao was obtained directly 
from plots of P/P^vs r/b; the second method in- 
volved measuring the slope of P/P«, plotted vs log 
(r/b), since 



It should be noted that much of tho P ft data gathers 
around the theoretical value of 0,715 for constant 
eddy viscosity. A comparison of these data with cal- 
culations of v^/p from the circulation profiles (see 
Eq, 17) indicated that the low values for P s were 


In the present experiments, quantitative axial 
velocity profiles were difficult to obtain, but tho 
visual studios revealed velocity dofects in the core 
at 4“ and 8° angles of attack. The most notable fea- 
ture of the axial flow was the presence of a lharp 
radial peak in the axial velocity profile near the 
vortex axis, This observation qualitatively agrees 
with the theoretical rosulls of Brown 8 regarding tho 
roll up of the momentum deficient boundary layer 
fluid into the vortex. Accordingly, the axial velo- 
city profile- whether representing a wake or a jet, 
assumes a Jyr dependence through the core. 

Measurements of the axial velocity on the axis 
were made by following fine lines of entrained hy- 
drogen bubbles which marked the vortex filaments in 
many cases. The magnitude of the maximum defect ve- 
locity depended strongly upon the angle of attack, 
as shown in Fig. 10, These measurements, furthermore, 
revealed that the velocity on the vortex axis accel- 
erated to a maximum before decaying with time, Accel- 
eration of the axial flow was accompanied by strong 
radial flow toward the axis in the core region. The 
mechanism for acceleration of the axial flow may be 
related to a stroamwise* pressure rise resulting 
from the decay of the integral for the radial pres- 
sure gradient: 



Since the stationary observer sees a velocity pro- 
file that decays with time, we have for the pressure 
gradient in fixed coordinates: 
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( 10 ) 



where we made the transformation z » U^t, The fluid 
surrounding the axis would ttas tend to accelerate 
toward she low pressure region behind the wing. With 
increasing time, the axial pressure gradient is di- 
minished, while viscous and turbulent friction cause 
the axial ilow to slow down, 


The Initial acceleration of the axial flow to- 
ward the wing at moderate angles of attack assumes 
significance in regard to the stability of the vor- 
tex core, that is, its resistance to turbulent eddies 
diffusing from the outer core region inward, betting 
the swirl velocity be invariant in the azimuthal di- 
rection and setting again z B U«t, the equation of 
continuity is given by 


1 3{nO I 3w 

7 ' Dr U. Di 
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Tf we assume a mean value for the axial acceleration 

in the coro 5 -y~> the radial inflow, u(r) , is approxi- 
mated by 


u(r) R* - 
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The momentum equation for the tangential flow is 
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Inside the core, so that by combining Eqs, 

(7) and (8) , We obtain for the rate of change of the 
tangential core velocity 

i (&.) (9) 

v ai u.d +w/uj \ at/ ^ 


z*i- u ill- 

Owen proposed a more sophisticated model, but in 
order to achieve agreement with experiment, also 
finally assumed that Vy is proportional to P M , Both 
models result in a rate of decay of v a proportional 
to t"i'‘. 

However, different rates of vortex decay have 
been previously observed, which, together with the 
data of this study, suggest that the eddy viscosity 
may vary with time, Particularly in towing experi- 
ments, where the airfoil is moved through an undis- 
turbed fluid, one may expect an initial increase of 
the eddy viscosity and a corresponding contribution 
to the rate of vortex decay, This reasoning seems 
logical in view of the fac- that the fluid is ini- 
tially free of turbulence, and that the rolled up 
boundary layer from the airfoil does not account for 
the turbulence levels found in the vortex field at 
some later time, In wind tunnel tests, and particu- 
larly in flight test situations where rsropulsien 
thrust contributes to the initial turbulence, one 
would expect a decay of the initial turbulence with 
time, and honce u decay rate for v a slower than 
t-1/2, billey's20 examination of the experimental 
flight data of nose and 0ee21 with respect to time 
varying eddy viscosity led to a decay of the maximum 
tangential velocity according to t"*' . In the work 
of Baldwin ot al,®on turbulent flow calculations, it 
was pointed out that the vortex initially is a non- 
equilibrium flow, i.e,, turbulence production f dis- 
sipation, Thus their calculated eddy viscosity de- 
creased sharply from the initial value £“ O.SxlO 3 ) 
and then reached a constant of Vj/i> <*> ISO, It is in- 
teresting to observe that during this time, the cir- 
cumferential velocity attained only a very small 
decay rate, As v,Jv approached a constant or pre- 
sumably an equilibrium value, the maximum tangential 
velocity reached the familiar dependence upon t "i/2. 


Hence, spin up of the core angular velocity is in- 
duced by acceleration of tbs axial core flow, The 
vortex filament is thereby strengthened and further 
stabilized. This phenomenon may be partly responsible 
for the reported persistence over long periods of 
time of well defined cores of aircraft trailing 
vortices, 

Effect of turbulent Viscosi ty on Vortex Decay 


The following describes an analytical solution 
for the decay of an originally potential vortex with 
turbulence when the eddy viscosity is allowed to 
vary with time, but remains constant in the radial 
direction. Since the actual developtaent of turbu- 
lence in the vortex flow is not considered, the so- 
lution shows the dependence of the rate of vortex 
decay upon particular assumed variations in n T (t) 
which are actually observed, 


Theoretical Considerations 

It is generally agreed that the effective, or 
turbulent eddy viscosity, present in the vortex 
is considerably higher than the molecular viscosity, 
and thus contributes to the rate of vortex decay, 
According to data from various investigators col- 
lected by Oven, estimates for v^/v fall between 
30 and 2000; depending upon the vortex Reynolds 
number, r^Jv. Kuhn and Nielsen, * 8 for instance, pro- 
posed a theory to calculate the vortex decay by 
using the z-coraponent of the eddy viscosity evaluated 
at the vortex axis. Their predictions fit the wind 
tunnel data of Chigier and Corsiglia 2 for values of 
the eddy viscosity ratio between 200 end 4SC, Squire^ 8 
suggested that Lamb's laminar solution can still be 
applied to tho turbulent vortex decay when the eddy 
viscosity is constant with time and the radial co- 
ordinate, and is given by 


We define the effective viscosity of the flow 

v t s v + T"* = T m j £H) 

whore is the initial eddy viscosity level in 

the flow after roll up. The symbol, T, is dimension- 
less time using appropriate parameters of the flow, 
e.g,, T " tfoo/c 2 , in order to solve the resulting 
differential equation, we make the further assump- 
tion that (Vj/v)# - ! 131 > 1, This assumption is justi- 
fied in wind tunnel and flight tests, except perhaps 
in low Reynolds number towing experiments. Under 
these circumstances, the two-dimensional vortex mo- 
tion is governed by 
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The solution to this equation is 



We note that the overall features of the flow 
remain the same as with constant v T ; hut the rate 
of change with tine now depends upon the develop- 
ment of the eddy viscosity, We want to establish 
the rate of change of the maximum circumferential 
velocity and the growth of the vortex core as m 
assumes various values, Solving for the radius at 
the point of maximum tangential velocity, a(t), and 
substituting r B a(t) in equation (13), we have 


of the vortex, tut had not reached the equilibrium 
level which is characteristic of the dependence upon 
This conclusion derives from the result that 
the rate of vortex decay remained proportional to 
t»7/8. por the rate of change of the eddy viscosity 
with tine, we obtain from Eq. (14) and tho data in 
Fig, S that m « 3/4, and thus (VrA')a ~ t3/4, Using 
the appropriate equations given in the papers of 
Baldwin et al. 9 and Donaldson* 2 we estimated the 
equilibrium eddy viscosity level to be (v^lv) e q U , » 
ISO, 

Experimental Results 

The variation of t » T /p with radius and time was 
computed from a graphical fit to the circulation 
profiles. Substituting 


v(a f t) 'w U. 



-!tn+ I)/l 


(14) 
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in Eq, (13) and solving for v^/v, we have 

MO „ (m + l)r s 

v " , . , r„ - rtr.t) 

K 4nt In “k — — — 


(17) 


In contrast to vortex decay with constant the 
growth of the vortex is proportional to tC^i)/ 2 , 
i.e,, time dependent changes in Vy contribute to the 
rate of vortex decay through the additional exponen- 
tial factor m/2. Table 1 demonstrates the signifi- 
cance of these results. 


Table 1, Time dependence of maximum circumferen- 
tial velocity and core radius upon 
assumed time variation of i>j/v. 
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Surveying the literature, including the results 
of this study, one finds that the range of time de- 
cay factors listed in Table 1 include those either 
previously measured or predicted using various 
theories of turbulent shear flow, The case m * -1 is 
of interest since turbulent vortex decay initially 
behaves remarkably in this way, namely with no ap- 
preciable decay of the circumferential velocity or 
growth of core radius, The initial value of m, of 
course, does not remain constant tut tends toward 
zero at the same rate at which equilibrium condi- 
tions are established in the vortex flow. On the 
other hand, the calculations of Baldwin et al. showed 
that the t"*/ 2 dependence occurred from the outset 
when the initial conditions were chosen correspond- 
ing to an equilibrium flow, The authors called these 
solutions self similar, but we recognize that all 
solutions with constant m possess the self similar 
property. Hence, m * 0 represents the law for turbu- 
lent vortex decay when equilibrium conditions lave 
been established in the flow. The available data on 
wake vortices suggest strongly that the time re- 
quired to reach this condition is proportional to 
the degree to which the initial flow conditions de- 
viated from those of equilibrium. 

In view of the findings above, our test results 
in a towing tank suggest primarily that the turbu- 
lence intensities were rising during the development 


Using m « 3/4 as established from tho data of Fig, 

5, profiles computed from Eq. (17) were accepted as 
representing levels of the eddy viscosity approxi- 
mately when Vj/v was invariant with r, or nearly so, 
at the location of the core radius. Surprisingly, 
this condition was satisfied in the majority of the 
examined vortex profiles, Some typical results are 
shown in Fig. 11, indicating an expected dependence 
upon Tbo through « arid IU , A summary of the eddy vis- 
cosity data obtained at the location of the core 
radips is plotted in Fig. 12, together with the wind 
tunnel data of Refs, 18, 22-25, It is seen that the 
eddy viscosity levels in the towing experiments were 
below those obtained in wind tunnel tests, 

Fig, 13 shows the tine dependence of eddy vis- 
cosity, (i'^/n) a , determined at the radius of maximum 
tangential velocity for the data taken at a * 8°, 

The indicated increase of {Py/v) a , proportional to 
t 3 / 4 , is in agreement with the valuo of m established 
from the slope in Fig, 5, 

Inte raction of the turbulent shear stress, r/p » 

- u'v', with the mean flow provides the significant 
mechanism for dissipation of the mean flow enorgy 
via production of turbulence, Nondlmensional turbu- 
lent shear stress profiles were computed from 

r/p JV _/9v v\ 

F. ! /b’ n?/b J \ar " 7 / (1«) 

Fig, 14 shows shear stress profiles for « « 8 s at 
various downstream distances. They agree with the 
computed results of Donaldsun and Baldwin et al, in 
shape and she relative location of the peak value 
with respect to the core radius. Rates of produc- 
tion of turbulent kinotic energy, normalized upon 
Foo /b 4 , were computed from 



Fig. IS shows representative profiles of turbulence 
production for the conditions of Fig, 14 in relation 
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to downstream distance, At the core radius, P attains 
approximately one-half its peak value, It is also 
evident that as the core radius grows, decreasing 
portions of the wean flow energy are dissipated In 
the core due to the increasing region of solid body 
rotation, 

£amaluiU»g..iifiaa£]is 

Using a towing tank facility, tangential veloc- 
ity profiles were measured for downstream of iv rec- 
tangular air foil. The resulting rates of decay of 
the maximum tangential velocity foAlowed a t* 7 '8 
dependence, and rates of growth of the vertex core 
were also greater than expected. Analysis of the 
data showed that this particular behavior was re- 
lated to the increase in the turbulent eddy viscosity 
with time; specifically, n T /n appeared to increase 
according to t 3 ' 4 . Values of circulation measured at 
the point of maximum circumferential velocity, al- 
though scattered, remained constant in the mean with 
downstream distance at about 0,67 Hw , High axial 
velocities with initial acceleration in the towing 
direction were observed and were attributed to a 
strongly stabilising influence upon the vortex 
filament. 

An analytical solution to the problem of vortex 
decay with time varying eddy viscosity revealed that 
only under equilibrium flow conditions does the self 
similar solution for the tangential velocity depend 
upon t'1/2. Wind tunnel and flight test data, as 
well as flow calculations for turbulent Vortex de- 
cay, have established these results as the expected 
long time behavior, Nonequilibrium flow with tem- 
porally decreasing turbulence intensity (eddy vis- 
cosity) yields other self similar solutions which 
decay slower tlian t _ l/ 2 . when Vy increases with time, 
as was tho case in the present experiments, the vor- 
tex decay proceeds at a rate faster than t "i/2, The 
low turbulence levels created by the airfoil in our 
tests resulted in conditions so far removed from 
equilibrium f.ow, that the eddy viscosity increased 
steadily wi :h time, thereby causing an accelerated 
rate of vortex decay. According to the data, the 
vortex could be considered decayed before it reached 
tho slower equilibrium decay rate, 
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Fig. 3 Schematic of hydrogen bubble technique in- 
dicating determination of tangential valoc- 
tty. 


Fig. 1 experimental set-up It. towing tank; 1-atrut, 
2-ailrror, 1 -window, 4-cameras, S-pulse gen- 
erator, 6-ground plate, 7-Hj bubbles, 8- 
alrfoll. 
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Fig. 2 Axial view of vortex field 
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rig. 7 Decay of circulation profiles with down- 
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Pig. 5 Decay of maximum tangential velocity. 
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Fig, S Similarity plot of circulation profiles, 
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Pig, 6 Growth of core radius behind airfoil 
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Fig, 9 Values of circulation determined at the 
care radius, 
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Fig, 10 Maximum axial velocity in towing direction 
after roll up. 
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Fig. 12 Comparison of experimental eddy viscosity 
data. 
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Fig, 11 Eddy viscosity profiles computed from ex 
perimental data and Eq. (17) . 


9 


Fig, 15 Time dependence of eddy viscosity deter- 
mined from data. 
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